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ABSTRACT 



The well-known self-similar solution for the two-point correlation function of 



CN ' the density field is valid only in an Einstein - de Sitter universe. We attempt 

^ to extend the solution for non - Einstein - de Sitter universes. For this purpose 

we introduce an idea of quasi-self-similar evolution; this approach is based on 
^ the assumption that the evolution of the two-point correlation is a succession of 

p i' stages of evolution, each of which spans a short enough period to be considered 

Q ■ approximately self-similar. In addition we assume that clustering is stable on 

■ scales where a physically motivated 'virialization condition' is satisfied. These 

■ assumptions lead to a definite prediction for the behavior of the two-point cor- 
relation function in the strongly nonlinear regime. We show that the prediction 

^ ■ agrees well with N-body simulations in non - Einstein - de Sitter cases, and 

■ discuss some remaining problems. 

Subject headings: cosmology: theory — large-scale structure of universe — grav- 
itation — methods: N-body simulations 



1. Introduction 

It is generally believed that the structure in our universe has grown out of tiny density 
fluctuations through gravitational instability. The growth of those density fluctuations is 
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completely described by linear theory on scales much larger than the correlation length 
of the density field (e.g., Peebles 1980). The behavior of the nonlinear density field on 
smaller scales, however, needs to be modeled with additional assumptions, and the resulting 
predictions should be verified and calibrated through extensive comparison with N-body 
simulations. 

The most popular prediction in the nonlinear regime is based on the combination of the 
self-similar evolution and the stable clustering ansatz (Peebles 1974, 1980; Davis & Peebles 
1977); if the linear power spectrum of density field follows a single power law oc A;", and the 
universe is described by the Einstein - dc Sitter model, the evolution proceeds in a self-similar 
manner since there is no characteristic scale in the system. This self- similarity, together with 
the ansatz that the clustering is stable in the strongly nonlinear regime, predicts that the 
logarithmic slope of the two-point correlation function ^ is equal to — 3(n -|- 3)/(n -|- 5). 

This solution has been widely applied in modeling the nonlinear gravitational cluster- 
ing (Hamilton et al. 1991; Peacock & Dodds 1994, 1996; Jain, Mo & White 1995) and in 
understanding the pair-wise velocity dispersions and thus the redshift-space distortion (Suto 
& Jing 1997; Jing, Mo & Borner 1998). In fact, the above prescription has been apphed 
even in the cases where the universe is not described by the Einstein - de Sitter model, 
and/or the linear power spectrum is not of a power-law form. The fitting formula for the 
nonlinear power spectrum by Peacock & Dodds (1996), for instance, takes account of the 
non power-law nature of the linear spectrum, but they have not examined the validity of the 
stable clustering solution in non Einstein - de Sitter models because 'if collapse occurs at 
high redshift, then Q = 1 may be assumed at that time' (Peacock & Dodds 1994). 

We revisit this issue in detail for the non - Einstein - de Sitter case, i.e., when the matter 
density parameter Qq is smaller than unity. Using N-body simulations with 64'^ particles, 
Suginohara ct al. (1991) already found that while the stable solution is reproduced well in 
the Einstein de - Sitter model, the slope of ^ becomes steeper in < 1 models (see also 
Suto 1993). This is not surprising at all because the derivation of the solution heavily relies 
on the scale-free nature of the Einstein - de Sitter model, that is, (1) a oc i^/^, where a is 
the cosmic scale factor, and (2) D{t) oc a, where D{t) is the hnear growth rate of density 
fluctuations. This motivates us to find a suitable modification of the self-similar solution 
which is also applicable to the non - Einstein - de Sitter case. 

In this paper we introduce an idea of quasi-self-similarity. This is based on the assump- 
tion that the evolution of the two-point correlation proceeds as a sequence of different quasi- 
self-similar stages, each of which is described by the locally self-similar solution determined 
by the cosmological parameters at that epoch. We also assume that the clustering becomes 
stable on scales where a 'virialization condition' is satisfied. These two assumptions lead to 
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a prediction of the behavior of the two-point correlation function in the strongly nonlinear 
regime. We show the extent to which the prediction agrees with the results of high-resolution 
N-body simulations. 



2. Quasi-Self-Similar Evolution 

2.1. The Self-Similar Solution in the Einstein-de Sitter Model 

For later convenience, we first briefly review the self-similar solution in the Einstein- 
de Sitter case (Peebles 1974, 1980; Davis & Peebles 1977). Suppose a system of particles 
with mass m. In the fluid limit, the system is described by the one-particle distribution 
function /(x, p,i), where x is the comoving coordinate, and p = ma^dx/dt is its canonical 
momentum. The distribution function /(x, p, t) obeys the Vlasov equation (Inagaki 1976; 
Peebles 1980): 

^^J!^^_m—— = (1) 

dt ma? dxi dxi dpi ' 

where the gravitational potential 4> satisfies 

VV = ^TiGma-^ [ fd^p. (2) 



In the Einstein - de Sitter model (a oc t'^^^), equation (1) admits a self-similar solution 
of the form: 

/(x, p, t) = r'^-'f{K/e, pA"+V3). (3) 

The value of the parameter a is determined by matching the solution in a linear regime as 
follows; first note that equation (3) implies the following form for the two-point correlation 
function ^{x, t): 

ax,t)^i(x/n. (4) 

If the initial power spectrum of density fiuctuations, Pinitiai(^), is proportional to A;", then 
^{x,t) at later epochs behaves as x~^"'~^^h'^^^ on large scales where linear theory is valid. 
Thus the value of a is exphcitly specified as 

4 

a — —. r. (5) 

3(n + 3) ^ ^ 
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In the small scale (nonlinear) limit, it is often assumed that the average proper separa- 
tion of pairs remains constant (the stable clustering hypothesis): 

{v2i{x,t)) = -ax, (6) 

where V21 is the relative peculiar velocity for a pair, and the brackets denote the average 
over pairs at a given comoving separation x. This assumption with equation (4) fixes the 
behavior of the correlation function in the nonlinear limit. Substituting equation (6) into 
the equation of the particle pair conservation: 

yields 

I = -A f [-'4 (8) 

ot ax^ox 

for ^ 3> 1. Thus ^ in the nonlinear regime should be of the form 

C^a^g{ax), (9) 

where g is an arbitrary function at this point. Finally the consistency with the above 
similarity solution requires that ^ should be given by the following power-law form: 

3(n+3) 4 3(n+3) 6 

^{x,t)ocx "+B i"+5 oc X "+5 a"+5. (10) 



2.2. Quasi-Self-Similar Clustering Solutions 

Now we attempt to generalize the above self-similar solution in a hypothetical universe 
where a oct^ and D{t) oc f^. In non-Einstein-de Sitter models, both p and q are not constant 
but change with time. As long as the time-scale of the change of p and q is smaller than 
that of the cosmic expansion, however, they can be regarded as constants for a short period 
at each epoch. 

We begin with the assumption of similarity for ^{x, t): 

ax,t) = ax/n. (11) 

Unlike in the Einstein - de Sitter case, any kind of self-similarity does not hold in a strict 
sense. However, the range of p and q considered here is close to that in the Einstein - de 
Sitter case, so it seems reasonable to assume that some form of self-similarity is realized 
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in an approximate sense. The simplest possibility is that, as in the previous subsection, 
^{x,t) has the form (11). With this ansatz we repeat the same procedure in the previous 
subsection, adopting Pinitiai(^) and the stable clustering hypothesis in the strongly 

nonlinear regime. 



In this case, one has ^ oc x ("+3)^2? linear regime, and equation (5) is replaced by 

2q 



a — 



n + 3' 



(12) 



Combining the form (9) for ^ in the strongly nonlinear regime yields, instead of equation 
(10), 



3(n+3) 



6/ 



3(k+3) 



^{x,t) oc X "+3+2/ a "+3+2/ = (ax) "+3+2/ a^, 



(13) 



for ,^ ^ 1. The quantity f = q/p in the above expression is in fact the familiar logarithmic 
derivative of D with respect to a. An excellent approximation to / is (Lahav et al. 1991) 



/ 



d\nD 
din a 



where 



and 



A(a) 



Q(a)- + M 



^0 



1 + 



n{a) 



Qo + {l-Qo- Ao)(a/ao) + Ao(a/ao)^ 

\o{a/aof 



+ (1 - ^0 - Ao)(a/ao) + Xo{a/ao)^ ' 
with Ao being the dimensionless cosmological constant at the present epoch oq. 



(14) 



(15) 



(16) 



2.3. Comparison with N-body Simulations 

Let us compare in detail the self-similar solution (eq. [10]) and the quasi-self-similar 
solution (eq. [13]) with high-resolution N-body simulations (Jing 1998; 2001b in preparation). 
The simulations consist of two Einstein - de Sitter models, two open models (Qq = 0.1, 
Ao = 0), and two spatially flat models (JIq = 0.1, Aq = 0.9). All the models employ scale- 
free initial power spectra Pinitiai(^) oc /c" with n = —1 and n — —2. Gravitational force 
calculation is based on the P^M algorithm. Simulation parameters are listed in Table 1, and 
further details of the simulations are described in Jing (1998, 2001b in preparation). 
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Figure la plots the two-point correlation functions which are appropriately scaled with 
respect to x according to the conventional self-similar solution (eq.[10]). It is clear that while 
Qq = 1 simulations are in good agreement with equation (10) beyond C, = 100 (indicated by 
an arrow for each model), = 0.1 models are not; the disagreement is particularly severe 
for the Ao = models. Figure lb shows the similar plot but with scaling for quasi-self-similar 
solution (eq.[13]), where we use the value of / evaluated at the present epoch. The figure 
indicates that the agreement with the simulations in Qq — 0.1 models indeed improves. 
Nevertheless the results in = 0.1 and Aq = 0.9 models do not yet show acceptable 
agreement with equation (13). The origin of the disagreement is studied by Suto, Taruya, 
& Suginohara (2001). 

This comparison points to the following two suggestions; (i) in < 1 models the 
conventional self-similar solution fails to describe the behavior of ^ for ^ > 100, and (ii) taking 
account of the dependence of the slope on the cosmological parameters as in equation (13) 
does not yet yield an acceptable prediction. With these points in mind we attempt to improve 
the quasi-self-similar solution in the next section. 



3.1. Virialization Condition and Scale Dependence of the Slope of ^ 

In the previous section we have applied the quasi-self-similarity only at the present 
epoch. In reality, however, the logarithmic slope of the correlation function: 



should be fixed locally at the epoch of the virialization of the corresponding scale, avir(a;). 
This naturally generates additional scale-dependence on the resulting solution through the 
time-dependence of /. More specifically, we attempt to incorporate this effect and to improve 
the model in the previous subsection as follows. 

In order to determine avir(a;) at each scale x, we need to make an assumption about the 
scale on which the system has been virialized at a given epoch a. Here we mainly adopt the 
following assumption: the system has just been virialized at a scale Xvir(a) where the volume 
average of the two-point correlation function reaches the critical overdensity of a virialized 
halo, Avir(a), predicted by the spherical collapse model: 



3. An Improved Model 



dine 



3(n + 3) 



(17) 



dlnx 



n + 3 + 2f{ayir{x)) 




(18) 
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where 

fy'^Wy- (19) 

^ Jo 

It seems natural to use ^ rather than ^ for the present purpose because the right-hand-side 
of equation (18) is the average density of a halo that has just been virialized. In Aq = 
models, the critical density Avir(a) is explicitly written as (Lacey & Cole 1993) 

cosh 77 = -^-1, (21) 
and its accurate fitting formula in Aq = 1 — f^o models (Nakamura & Suto 1997) is given by 

r r I -| 0.9052 ^ 

Avir(a) ~ IStt^I 1 + 0.4093 -ij |. (22) 

The 'virialization condition' (18) may be somewhat arbitrary, but is perhaps most natural 
and physically motivated among other choices. 

Furthermore we assume that, for x < Xyir{a), the stable clustering condition is satisfied, 
i.e., ^(x, a) is of the form (9). This is equivalent to saying that, at a fixed physical scale 
r = ax, the slope of ^ is kept constant while its amplitude grows in proportion to a^ (c.f., 
eq.[13]). 

The above assumptions are not yet sufficient in predicting ^{x) for a given model. The 
amplitude of ^ at an arbitrary point needs to be specified by hand because our model does 
not predict the overall amphtude of ^{x). Once the value of ^ is fixed at a scale x — e 
(for example, in the next subsection, we take ^grav, the gravitational softening length, as 
e and use the value of ^(Sgrav) in the simulations), we can compute ^{x) at a — ao up to 
Xvir,o = Xyir{ao). The procedure is as follows: 

1. Set boundary condition. — Below an extremely small scale xi, where ^ is sufficiently 
large, our model prediction should reduce to the conventional self-similar solution in the 
Einstein - de Sitter universe; i.e., 

C OC ^-3("+3)/(n+5) foj. ^ < (23) 

This is because the epoch of the virialization of the corresponding scale is sufficiently early, 
the universe is indistinguishable from the Einstein - de Sitter model (Peacock & Dodds 
1994). Equation (23) implies that 

e>0 = ^-^^{xi). (24) 
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We first choose xi arbitrarily and set the value of ^{xi) to be sufficiently large, say, 10^. This 
is the starting point of our computation. 

2. Compute ayir{x). — Prom the above assumptions it follows that 



At the current scale x, we solve equation (25) for avir(a;). 

3. Advance x . — Substituting the value of a^u-ix) into equation (17) gives the local slope 
d\n^/d\nx at x. Then we can advance a; by a small interval Ax and compute ^{x + Ax), 
and then ^{x + Ax) . 

4. We repeat the procedures 2 and 3 until ^{x) becomes equal to Avir(ao), which 
corresponds to the present virialization scale Xyir.o- 

5. Normalization. — Finally we shift the resulting solution so that it can match the 
given amplitude at x = e, using the fact that, if ^'■^^(x) is a solution to equation (17), so is 
^(^^(x) = ^^^\ax) with an arbitrary constant a. 



Figure 2 compares our improved model predictions with the N-body results. Using 
the scaling relation described above, we match the amplitude of our solution to that of 
the simulations at a; = Sgrav for each model (the innermost symbols). The length scale x is 
normalized by a;vir,o = a;vir(ao) and we show the results only in the virialized regime, x < Xvir.o- 
Clearly, our predictions are in good agreement with simulations for ^ > 200 in all models. 
In particular the predicted dependence on the spectral index n is excellently reproduced in 
the simulation results. Given the simphcity of our procedure, this may be regarded as a 
considerable success. For ^ < 200, however, the slopes of our predicted correlation functions 
are shallower than those of simulations (Figs. 2a and 2b). 

Let us compare our model predictions with the conventional self-similar solution (eq. [10]) 
in more detail. The power-law slope of the conventional self-similar solution is shown by the 
dot-dashed lines in Figure 2. In Qq = 0.1 and Aq = 0.9 cases both our model and the 
simulation have roughly the same slope as the conventional one near x — £grav) but in both 
of them the slope becomes steeper as the scale approaches to a;vir,o- More impressively, in 
Qq = 0.1 and Aq = 0.0 cases both our model and the simulation have a steeper slope on all 
scales shown in the figure than the conventional one. 




(25) 



3.2. Comparison with N-body Simulations 



Also shown in Figure 2 are the correlation functions obtained by Fourier-transforming 
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the fitting formulae for the nonlinear power spectra by Peacock & Dodds (1996). The 
Peacock-Dodds formula is extremely useful because it gives not only the shape but also the 
amplitude of the two-point correlation function from linear to strongly nonlinear regimes. 
The agreement between the Pcacock-Dodds formula and the simulations is very well except 
for f^o = 0.1 and Aq = models. In these cases the Peacock-Dodds formula systematically 
underpredicts the simulation results. Nevertheless this could be adjusted somehow by shifting 
the amplitude of the Peacock-Dodds formula so as to match the simulation at x — £grav 
(see the dotted hues in Fig. 2). Rather an important advantage of our model over the 
Peacock-Dodds prescription is that it does successfully predict the slope of ^ up to a scale 
where deviation from the conventional one, — 3(n -|- 3) / (n -|- 5) , is significant. Note that the 
slope of ^ computed from the Pcacock-Dodds formula is, on scales where it deviates from 
— 3(n + 3)/{n -\- 5), simply an interpolation from the numerical simulations. In this sense, 
our result implies that there is room for further improvement in the original Peacock-Dodds 
formula, and our present model may be useful for that purpose. 

One may wonder if it is possible to improve our model predictions by varying the 
condition (18) somehow. Figures 3 and 4 present those results. The top panels in those 
figures adopt ^{x^iT.{a)) = 2A^i^{a) instead of equation (18). While the agreement between 
the quasi-self-similar prediction and the simulations is indeed improved for = 1 models, 
this is not the case for Qq < 1 models. 

The middle and bottom panels show the results for ^{x^ij.{a)) = 100, and ^{x^i^{a)) = 
300, respectively, instead of equation (18). In each case we obtain acceptable agreement for 
both flo = 1 and Qq < 1 models simultaneously, although the modified conditions lose the 
physical basis and should be regarded as empirical at best. 



3.3. Validity of the stable clustering hypothesis 

The stable clustering hypothesis is an essential ingredient in our model, but has some- 
what been in doubt in the recent literature (Yano & Gouda 2000; Ma & Fry 2000; Caldwell 
et al. 2001). We argue here, nevertheless, that the hypothesis still remains to be a reasonable 
assumption. 

Yano & Gouda (2000) relate the inner density profile of virialized halos with the velocity 
parameter h = —{v2i)/(ix. They claim that h should approach zero in the nonlinear limit if 
the logarithmic slope of the inner density profile is larger than —3/2. As pointed out by Ma 
& Fry (2000), their argument is valid only when all halos have an equal mass, and should 
be completely altered for a realistic mass function. 
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On the basis of simulations by Jain (1997), Caldwell et al. (2001) find for a variety of 
cosmological models a universal relation between and h, and propose a fitting formula 
describing the relation. Although extrapolating this formula implies that /i — in the 
nonlinear limit, they claim that their formula is valid for ^ 10^- Thus their results are 
not inconsistent with the idea that /i ~ 1 in the strongly nonlinear regime. In fact, more 
recent numerical work supports the stable clustering assumption (Jing 2001a; Fukushige & 
Suto 2001). 

4. Discussion 

We have shown that the conventional self-similar solution in the Einstein - de Sitter 
universe does not describe the behavior of two-point correlation functions in l^o 7^ 1 models 
for strongly nonlinear regimes of cosmological interest, ^ < 10"^. Instead, we have proposed 
a simple model to describe the two-point correlation functions in strongly nonlinear regimes 
by introducing the quasi-self-similar ansatz. In fact we have shown that the resulting model 
predictions in non - Einstein - de Sitter universes agree better with the high-resolution 
N-body simulations. 

On the other hand, our current model is not fully successful yet in the sense that the 
predicted behavior for ^ < 200 systematically differs from that observed in simulations. 
Empirically this situation can be improved by an appropriate choice of the virialization 
threshold ^vir- While the physical meaning of this procedure is not clear, this may be related 
to some other physics that we omit in the present simple prescription. After all the regime 
for ^ < 200 may not be completely dominated by the stable clustering evolution, and the 
effect of linear and quasi-linear evolution is likely to be important as well. In fact it may be 
the case that the stable condition (6) is not realized instantaneously when the virial condition 
(18) is satisfied, and that the slope at the corresponding scale approaches the value implied 
by equation (17) only gradually as the scale becomes more strongly nonlinear (for instance, 
^ > 200) and completely decoupled from the linear evolution of the environment. We plan 
to check this interpretation using the time-evolution of the simulation results in due course. 

One of the most important applications of the conventional self-similar solution is the 
fitting formula for the nonlinear power spectrum by Peacock & Dodds (1996), which is 
based on an idea originally proposed by Hamilton et al. (1991). In the strongly nonlinear 
regime, on which we focus in this paper, the Peacock-Dodds formula agrees well with the 
simulations for Oq = Ij and also for fio < 1 and f^o -|- Aq = 1. We have found, however, that 
in fio = 0.1 and Aq = cases the formula fails to fit the simulations. This may be understood 
as follows. In such cases the slope of ^ asymptotically approaches that in conventional self- 
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similar solution only on scales where ^ is extremely large. This means that the gap between 
these asymptotic scales and linear scales is rather big, and it cannot be simply interpolated. 
In contrast our model predicts the shape of ^ up to the scale corresponding to ^ ~ a few 
hundred, so interpolation between this scale and the linear scale may be much easier. Thus 
our model may be useful in improving the Peacock-Dodds formula, especially in Qq < 1 and 
Ao = cases. 

In summary, although our proposed model still needs to be improved, the degree of 
success for ^ > 200 is encouraging, and useful at least empirically. We hope that the idea 
of quasi-self-similar evolution may give a useful insight towards a better understanding of 
nonlinear dynamics of the mass density field in the universe. In particular, we plan to 
improve the existing fitting formula for the nonlinear power spectrum (Peacock & Dodds 
1994, 1996) by applying the quasi-self-similar ansatz. 
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acknowledges support from Research Fellowships of the Japan Society for the Promotion 
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Education, Science, Sports and Culture of Japan (07CE2002, 12640231) to RESCEU, and 
by the Supercomputer Project (No.99-52, No.00-63) of KEK. 



- 12 - 
REFERENCES 

Caldwell, R. R., Juszkiewicz, R., Steinhardt, P. J., & Bouchet, F. R. 2001, ApJ, 547, L93 
Davis, M. & Peebles, P. J. E. 1977, ApJS, 34, 425 
Pukushige, T. & Suto, Y. 2001, ApJ, 557, Lll 

Hamilton, A. J. S., Kumar, P., Lu, E. & Matthews, A. 1991, ApJ, 374, LI 
Inagaki, S. 1976, PASJ, 28, 77 
Jain, B. 1997, MNRAS, 287, 687 

Jain, B., Mo, H. J. & White, S. D. M. 1995, MNRAS, 276, L25 

Jing, Y. P. 1998, ApJ, 503, L9 

Jing, Y. P. 2001a, ApJ, 550, L125 

Jing, Y. P., Mo, H. J. & Bdrncr, G. 1998, ApJ, 494, 1 

Lacey, C. & Cole, S. 1993, MNRAS, 262, 627 

Lahav, O., Rees, M. J., Lilje, P. B. & Primack, J. R. 1991, MNRAS, 251, 128 

Ma, C.-P. & Fry, J. N. 2000, ApJ, 538, L107 

Nakamura, T. T. & Suto, Y. 1997, Prog. Theor. Phys., 97, 49 

Peacock, J. A. & Dodds, S. J. 1994, MNRAS, 267, 1020 

Peacock, J. A. & Dodds, S. J. 1996, MNRAS, 280, L19 

Peebles, P. J. E. 1974, ApJ, 189, L51 

Peebles, P. J. E. 1980, The Large-Scale Structure of the Universe (Princeton: Princeton 
Univ. Press) 

Suginohara, T., Suto, Y., Bouchet, F. R. & Hernquist, L. 1991, ApJS, 75, 631 
Suto, Y. 1993, Prog. Theor. Phys., 90, 1173 
Suto, Y. & Jing, Y. P. 1997, ApJS, 110, 167 

Suto, Y., Taruya, A. & Suginohara, T. 2001, Progress in Astronomy (Supplement), 19, 67 



- 13- 

Yano, T. & Gouda, N. 2000, ApJ, 539, 493 



This preprint was prepared with the AAS IM^jX macros v5.0. 



-14- 



Tablc 1. Simulation parameters. 
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Fig. 1. — Two-point correlation functions in the simulations scaled according to (a) conven- 
tional self-similar solution, and (b) quasi-self-similar solution with scale independent slope. 
The solid lines arc the results for fio = 1 models, and the symbols are those for Qq = 0.1 
models with either Aq = or Aq = 0.9. Error bars are estimated from three realizations 
for each model; only the error bars for n = —2 and (i^o? -^o) = (0.1,0.9) are shown. Other 
models have smaller error bars, which are not shown for clarity. The arrows correspond to 
^{x) = 100 for each model. 
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Fig. 2. — Quasi-self-similar solutions with a virialization condition ^(xvir(a)) = Avir(ci). The 
solid lines are obtained by numerically integrating eq. [17]. The length scale x is normalized 
by iCvir 0- The dashed lines are computed using the Peacock-Dodds formula. The dotted 
lines also show the Peacock-Dodds formula, but this time each result has been horizontally 
shifted to match the simulation ai x — £grav The dot-dashed hnes show the slope of the 
conventional self-similar solution. The symbols are the results of the simulations, (a) n = — 1 
and (Jlo,Ao) = (1.0,0.0). (b) n = -2 and (llo,Ao) = (1.0,0.0). (c) n = -1 and (Jlo,Ao) = 
(0.1,0.9). (d) n = -2 and (Qo,Ao) = (0.1,0.9). (e) n = -1 and (QoAo) = (0.1,0.0). (f) 
n = -2 and (Qq, Aq) = (0.1, 0.0). 
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Fig. 3. — Quasi-self-similar solutions forn = —1. The solid lines are obtained by numerically 
integrating eq. [17]. The symbols are the results of the simulations, (a) Qq = 1 model with 
a virialization condition ^(a;vir(a)) = 2Avir(a). (b) < 1 models with ^(a;vir(a)) = 2Avir(o). 
(c) fio = 1 model with ,^(a:vir(a)) = 100. (d) < 1 models with ,^(,Xvir(a)) = 100. (e) l^o = 1 
model with ^(a;vir(a)) = 300. (f) flo < 1 models with ^{xyij:{a)) = 300. 
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Fig. 4. — As for Fig. 3, but for n — 



-2. 



